Nuclear matrix elements for double-/? decay 
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Background: Direct determination of the neutrino mass through double-/? decay is at the present 
time one of the most important areas of experimental and theoretical research in nuclear and particle 
physics. 

Purpose: We calculate nuclear matrix elements for the extraction of the average neutrino mass in 
neutrinoless double-/3 decay. Methods: The microscopic interacting boson model (IBM-2) is used. 
Results: Nuclear matrix elements in the closure approximation are calculated for ''*^Ca, ^^Ge, ^'^Se, 



^Zr, 



'Mo, 



'Pd, 



^-Cd, i^^Sn, 



'Te, 



"Te, 



^Nd, 



"Nd, ''^"Sm, ^""Gd, and ^-^^Pt decay. 



Conclusions: Realistic predictions for the expected half-lives in neutrinoless double-/? decay with 
light and heavy neutrino exchange in terms of neutrino masses are made and limits are set from 
current experiments. 



PACS numbers: 23.40.Hc,21.60.Fw,27.50.+e,27.60.+j 



I. INTRODUCTION 

In recent years, the possibility of a direct measurement 
of the average neutrino mass in neutrinoless double-/? de- 
cay has attracted considerable attention. Three scenarios 
have been considered [H-HI, shown in Fig.[T] After the dis- 
covery of neutrino oscillations 0^3, attention has been 
focused on the first scenario (a). In very recent years, 
the second scenario (b) has again attracted attention Q. 
For all three processes (Oj//3/3, QvhPfi, and Oi//3/3M), the 
half-life can be factorized as 



(1) 



where Goi/ is a phase-space factor, Mojy is the nuclear 
matrix element, and f(rni,Uei) contains physics beyond 
the standard model through the masses nii and mixing 
matrix elements Uei of neutrino species. 






Figure 1. Neutrinoless double-/? decay mechanism for (a) light 
neutrino exchange, (b) heavy neutrino exchange, and (c) Ma- 
joron emission. 

In addition to the neutrinoless modes, there is also the 
process allowed by the standard model, 2i//3/3, depicted 
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in Fig. [5] For this process, the half-life can be, to a good 
approximation, factorized in the form 



'1/2 



(2) 



(The factorization here is not exact and conditions under 
which it can be done are discussed in Ref. Q and Sec. 
III). 

A-2 




n A-2 



Figure 2. Double-/? decay mechanism with the emission of 2v. 

The processes depicted in Figs. [1] and [2] are of the type 
{A, Z) ^ {A,Z + 2) + 2er + anything. (3) 
In very recent years, interest in the processes 



{A, Z) {A,Z -2) + 2e+ + anything 



(4) 



has also arisen. In this case there are also the com- 
peting modes in which either one or two electrons are 
captured from the electron cloud {Ov^EC, 2^/3EC, and 
2j/ECEC). Also for these modes, the half-life can be fac- 
torized (either exactly or approximately) into the prod- 
uct of a phase-space factor and a nuclear matrix element, 
which then are the crucial ingredients of any double- /3 de- 
cay calculation. 

Recently, we have initiated a program for the system- 
atic evaluation of both quantities. The evaluation of 
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the phase-space factors (PSFs) for Oi//? (3 and 2i/^ (3 
has been reported in [H and that for QvP'^P'^ , 2I//3+/3"'", 
Oi^/3+EC, 2J//3+EC, and 2z/ECEC is in preparation 0. 
The main difference betveen this new calculation of PSFs 
and older standard approximations is a few percent for 
light nuclei {Z = 20), about 30% for Nd [Z = 60), and a 
rather large 90% for \J {Z = 92), the correction increas- 
ing as a power of Za. In this article, we concentrate on 
nuclear matrix elements Mq^, and M21, for Oi//3~(3~ and 
2vj3~ I3~ . Calculations of the nuclear matrix elements for 
positron emission have also been completed and will be 
reported in a subsequent paper fioj . Nuclear matrix ele- 
ments have been evaluated in a variety of models, most 
notably the quasiparticle random phase approximation 
(QRPA) and the interacting shell model (ISM). Results 
up to 1998 are reviewed in Refs. and (T^. In 1999 
a new formulation of Qvfip was introduced j3| and cal- 
culations within the QRPA jUl and the ISM Q were 
performed, as well as within other models, as discussed 
in the following Sec. II C. In 2009, we developed [HI 
a new method to evaluate nuclear matrix elements for 
double-/? decay within the framework of the microscopic 
interacting boson model (IBM-2). The advantage of this 
method is that it can be used in any nucleus and thus all 
nuclei of interest in both 13^/3^ and /?^/?^ decay can be 
calculated within the same model. 

The calculation of the nuclear matrix elements is done 
in the closure approximation. This approximation is 
good for Quj3p decay, since the average neutrino momen- 
tum is of the order of 100 MeV/c. It is, in principle, not 
good for 2i//3/3, since the average neutrino momentum is 
of the order of few MeV/c. However, formally the ap- 
proximation is still valid if one appropriately defines the 
closure energy. The advantage of the closure approxima- 
tion is that all calculations for the processes depicted in 
Figs. [1] and [2] can be done simultaneously, by changing 
the so-called neutrino potential, as discussed in the sec- 
tions below, thus eliminating systematic (and accidental) 
errors in the calculation, especially in the ratio of matrix 
elements for different processes. 



In this article, we report the results of our calculations 
for the nuclei listed in Table HI A selected number of de- 
cays wer considered in (isj and preliminary results were 
presented in (Tg], [T3| • Here we report the complete list of 
results divided into Ovfi/S (light neutrino exchange) and 
Qvhl3f3 (heavy neutrino exchange). Sec. II, and 2i^/3/3, 
Sec. III. By using these results we also set some limits 
on the mass of light {m,y) and heavy (m^^) neutrinos. 

II. NEUTRINOLESS DOUBLE/3 DECAY (Ou/SfS) 

A. Transition operator 

The theory of Oj^/3/3 decay was first formulated by 
Furry [l8| and further developed by Primakoff and Rosen 



Table I. Double-/3 decays considered in this article, their Q- 
values, and their isotopic abundances. 



P P transition 


Qflfl(keV) 


P(%) 


2oCa28 - 


->-22Ti26 


4272.26 ± 4.04 


0.187 ± 0.021 


32Ge44 - 


■>34Se42 


2039.061 ± 0.007 


7.73 ± 0.12 


34Se48 — 


^36Kr46 


2995.12 ± 2.01 


8.73 ± 0.22 


4oZr56 — 


^!iM054 


3350.37 ± 2.89 


2.80 ± 0.09 


i§°M058 


— ^44^R-U56 


3034.40 ± 0.17 


9.82 ± 0.31 


46''Pd64 


-^■48"Cd62 


2017.85 ± 0.64 


11.72 ± 0.09 


ig®Cd68 


->-5o'^Sn66 


2813.50 ±0.13 


7.49 ±0.18 


124o 

50 Bn74 - 


->'52 1672 


2286.97 ± 1.53 


5.79 ± 0.05 


52 -1-676 - 


->-54 Ae74 


865.87 ± 1.31 


31.74 ±0.08 


130rp„ 

52 -1-678 - 


->-54°Xe76 


2526.97 ±0.23 


34.08 ± 0.62 


54 Ae82 ■ 


-^■ie^Bago 


2457.83 ± 0.37 


8.8573 ± 0.0044 


I^Ndss 


— >-62*Sm86 


1928.75 ± 1.92 


5.756 ±0.021 


60 Ndgo 


— >-62''Sm88 


3371.38 ± 0.20 


5.638 ± 0.028 


62''Smg2 


— ^M^'Gdgo 


1251.03 ± 1.25 


22.75 ± 0.29 


64°Gd96 


, i60n,^ 

^■66 ^794 


1729.69 ± 1.26 


21.86 ±0.19 


198 pi 
78 -rtl20 


, 198tt„ 
^80 ngii8 


1047.17 ± 3.11 


7.36 ±0.13 



[l9| , Molina and Pascual |20| , Doi et al. [ij , Haxton and 
Stephenson [2l|, and, more recently, by Tomoda Q and 
Simkovic et al. Q. All these formulations often differ 
by factors of 2, by the number of terms retained in the 
nonrelativistic expansion of the current and by their con- 
tribution. In order to have a standard set of calculations 
to be compared with the QRPA and the ISM, we adopt 
in this article the formulation of Simkovic et al. Q . The 
transition operator in momentum space, p = is writ- 
ten as 

T[p)^H[p)f{m,,U,,) (5) 
where for light neutrino exchange 

/(m„C/eO = ^, M= E {U^ufmu, (6) 

k— light 

while for heavy neutrino exchange 

k— heavy 

The (two-body) operator H{p) can be written as 
Hip) = J2 ^nrl [-h^iP) + h^''{p)<^n ■ aw 

n,n' (8) 

with the tensor operator defined as 

^nn' ^ 3 [((T„ • p) {an' ' p)] - (?„ • CT„' . (9) 

The Fermi (F), Gamow-Teller (GT), and tensor (T) con- 
tributions are given by 

h^{p) = h^yip) 

h'^^ip) = h'ilip) + h%{p) + h%{p) + h'iJM (10) 
h^{p) = f^Apip) + hi pip) + hlj^fip). 
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The terms AP. PP, and MM are higher order corrections 
(HOC) arising from weak magnetism (M) and induced 
pseudoscalar terms (P) in the weak nucleon current. The 
terms h^-^'^-'^{p) can be further factorized as 



F,GT.T 



(p) =v{p)h 



F.GT.T 



(P) 



(11) 



where v{p) is called the neutrino potential and are the 
f^F,GT,T^p^ the form factors. A list of form factors is 
given in Ref. Q and recast in the form used by us in 
Table |IT1 In this table, the finite nucleon size (FNS) is 
taken into account by taking the coupling constants gv 
and QA as momentum dependent 



ffv(p^) = 9v 
QAip^) = 9 A 



1 + ^ 



(12) 



1 + -K 



The value of My is well fixed by the electromagnetic form 
factor of the nucleon, = 0.71(GeV/c^)2 Haj and gv = 
1 by the hypothesis of conserved vector current (CVC). 
The value of Ma is estimated to be Ma = 1.09(GeV/c2) 
[H and gA = 1.269 



Table II. Form factors in the formulation of Q adapted to 
our calculation, rrip and rrin are, respectively, the proton and 
pion mass and = 3.70 is the isovector anomalous magnetic 
moment of the nucleon. 

Term h{p) 



2 (at./ai) 
9 a 



tGT 



/■AP 



(l+pVAfi)" 
3 (i+pVW^)" 



^2 [2.?. 



T 

■AP 
T 

PP 



2 2 



2 9v 



— "-AP 
TGT 

— llpp 

17 GT 
2 "A/i\/ 



The neutrino potential v(p) is given, in the closure ap- 
proximation, for light neutrino exchange by 



v(p) = -- 



'^p(p + Aj 



(13) 



For heavy neutrino exchange, the neutrino potential is 
given by 



v{p) 



(14) 



The contributions in momentum space, h^''^'^''^ (p), 
can be converted to the contributions in coordinate space, 
l-^F,GT,T by taking the Fourier-Bessel transforms 



F,GT,Ti 



(r) =- / jxipr) 



P{P + A] (15) 



X p^dp, 

for light-neutrino exchange and 



FnT,T 



(r) 



2 

1" ./o 

„2 



j,(pr)^F.«^.^(p) 



(16) 



X p dp, 



for heavy neutrino exchange. Here A = for Fermi and 
Gamow- Teller contributions and A = 2 for tensor contri- 
butions. 

Finally, an additional improvement is the introduction 
of short-range correlations (SRC). These can be taken 
into account by multiplying the potential V{r) in coor- 
dinate space by a correlation function /(r) squared. The 
most commonly used correlation function is the Jastrow 
function 



fj{ 



1 



(17) 



with a = 1.1 fm^^, b = 0.68 fm^^ and c = 1 for the 
phenomenological Miller-Spencer parametrization (25j, 
and, in recent years, the Argonne/CD-Bonn parametriza- 
tions (Hi a = 1.59/1.52 fm~^ b = 1.45/1.88 fm"^ and 
c = 0.92/0.46. Since our formulation is in momentum 
space, we take into account SRC by using the Fourier- 
Bessel transform of fj{r). 

In assessing the "goodness" of Simkovic's formulation 
it is of interest to compare it with Tomoda's formulation. 
Apart from some differences in definitions, namely the 
fact that Tomoda defines the transition operator with a 
factor of 1/2 in front of Eq. (8), see Eq. (3.31) of Ref. 0, 
and the tensor operator with a factor of 1/3 in front of 
Eq. (9), see Eq. (3.54) of Ref. Q and with a plus sign in 
front of the tensor operator in Eq. (8), in contrast with 
Eq. (13) of Ref. [I, and a nuclear radius R = tqA^^^ with 
ro = 1.2 fm instead of 7'o = 1.1 fm of differences which 
have caused, however, considerable confusion in the lit- 
erature, the main difference between Tomoda's formula- 
tion and Simkovic's formulation is that Tomoda considers 
more terms in H{p), nine in all, three GT terms, three F 
terms, one T term, one pseudoscalar (P) term, and one 
recoil (R) term. Also, except for the terms hvv{p) and 
hAAip), where the form factors and potentials coincide, 
all other form factors and potentials in Tomoda's formu- 
lation are different from those in Simkovic's formulation. 
Although in this article we report results in the latter for- 
mulation, we note that we have results available for seven 
of the nine terms in Tomoda's formulation, the three GT 
terms, the trhree F terms, and the T term. The form 
factors and potentials for these seven terms are listed in 
Table VIII of our Ref. jl5|. In Table II of the same ref- 
erence, we also show that the contribution of additional 
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terms Xqt'jXf' ^^'^ x't ^'^^ negligible and thus the 
assessment of the "goodness" of Sinikovic's formulation 
must reflect this point. 



Matrix elements 



We consider the decay of a nucleus into a nucleus 
2+2^^-2- An example is shown in Fig. [3l The nuclear 
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Figure 3. The decay 32Ge44 
decay. 



>34Se42, an example of double-/? 



matrix elements are those of the operator H{p) of Eq. (8) 
between the ground state of the initial nucleus and the 
final state with angular momentum Jp 



Mo, = (^^X;Ot\H{p)\^Y;JF 



(18) 



If the decay proceeds through an s-wave, with two leptons 
in the final state we cannot form an angular momentum 
greater than one. We therefore calculate, in this article, 
only Oi/P/S matrix elements to final 0^" states, the groimd 
state of, and the first excited state 0^, for which in a 
previous article Q we have calculated the phase-space 
factors. The form factors in Table |IT] have a common 
factor of g^, except hyy. They depend on and gy. 
We write 



9a 



2 j^(Ou) 



— I M 

9A 



(0^) 



T ' 



(19) 



with the ratio gy / g^ explicitly displayed in front of 
The ratio gv/gA is also implicitly contained in 



M 



F 

AOu) 



Mq^ and Mj. " through the terms h^ff^,j and /ima/ (see 



Table |TT1) , and the matrix elements Mq^^-* , Mp 



and 



are defined as 



M, 



GT 
M 



(^X:0+|/^«^(p)/gi|V:J^), 
f^^^{^X;Ot\h^{p)/g'y\^Y;JF), 



(20) 



{Oi'} _ A 



X-OXW{v)lg\\ Y;Jf 



The reason for this separation is that the calculated 
single-/? decay matrix elements of the GT operator in 
a particular nuclear model appear to be systematically 
larger than those derived from the measured ft values 
of the allowed GT transitions. The simplest way of tak- 
ing into account this result is by introducing an effective 
gA.eff, also sometimes written as gA,eff = 95a, where 
g is a quenching factor. The quenching of gA will be 
discussed in Sec. IIIII Here we report results of the calcu- 
lation of M^O") in IBM-2 with the free values = 1 and 
gA = 1.269. These form the baseline for any discussion 
of the nuclear matrix elements (NME) in Qvfip. 

In order to evaluate the matrix elements we make use 
of the microscopic interacting boson model (IBM-2) 27 . 
The method of evaluation is discussed in detail in 15 . 



Here we briefly mention the logic of the method, which 
is a mapping of the fermion operator H onto a boson 
space and its evaluation with bosonic wave functions. 
The mapping [2^ can be done to leading order (LO), 
next to leading order (NLO), etc. In Ref. [l^ we showed, 
by explicit calculations, that NLO terms give, in gen- 
eral, negligible contribution, < 1%. In this article, we 
present only LO calculations. The matrix elements of 
the mapped operators are then evaluated with realistic 
wave functions, taken either from the literature, when 
available, or obtained from a flt to the observed ener- 
gies and other properties {B{E2) values, quadrupole mo- 
ments, B{M\) values, magnetic moments, etc.). The val- 
ues of the parameters used in the calculation are given 
in the appendices. In Appendix A, we give the neutrino 
potential and its parameters. In Appendix B, we list the 
single-particle and -hole energies and strengths of inter- 
action. In Appendix C, we give the parameters of the 
IBM-2 Hamiltonian for each nucleus considered in this 
article, together with their references. As shown in the 
references quoted in Appendix C, the quality of the IBM- 
2 wave functions ranges from very good to excellent for 
nuclei with ^ > 70 where collective features are very 
pronounced, especially in deformed nuclei. As an exam- 
ple, we show in Fig. S] a comparison between calculated 
and experimental spectra for the pair of nuclei ^^°Nd and 
I'^Sm. For nuclei with A < 70, the IBM-2 description is 
only approximate, and one needs to go to the isospin con- 
serving versions IBM-3 jlgj and IBM-4 j30|. Nonetheless, 
we will report, for the sake of completeness, also results 
for ''^^Ca decay, with the proviso that these are rather ap- 
proximate. Also, in some cases, intruder conflgurations 
play a role, especially in the structure of the excited 0"*" 
state, and one needs to go to the conflguration-mixing 
version IBM2-CM [SlJ. All these improvements will be 
reported in subsequent papers. 
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1T27 



2+ / 


i086 


2* 


848 


6+ 


708 


0+ / — 


\^ 


4+ 


374 


2+ 


134 


0+ 






"Nd 



4* 


1138 


8^ 


/l30 




1062 


1* 


851 


6+ 


720 


0+ / — 


\675 


4+ 


381 


2+ 


130 


0+ 






"Sm 



4+ 


1614 


T 


1423 


6+ 


1272 


T 


1049 


0+ 


813 


4^^ — 


^740 


2+ 


314 


0+ 






4+ 


1449 


6+ 


1279 


2^^ 


1194 


2+ 


1046 


4+ 


773 


— 


/740 


1* 


334 


0+ 






exp 



exp 



Figure 4. Comparison between calculated and experimental low-lying spectra for the pair of nuclei Nd and Sm. 



C. Results 

The matrix elements of the operator have di- 

mension fm~^. It has become customary to quote the 
values of M^""^ by multiplying by the nuclear radius in 
fm, R = RqA^I^ , with i?o = 1-2 fm. The matrix elements 
are then dimensionless. 



Ov^ifi decay with light neutrino exchange 



In Table IIII| we show the results of our calculation 
of the matrix elements to the ground state, 0^, and 
first excited state, Oj, broken down into GT, F, and 
T contributions and their sum according to Eq. (19). 
We note that since we are covering all nuclei from 
A ~ 76 to ^ = 198, we have two classes of nu- 
clei, those in which protons and neutrons occupy the 
same major shell {A = 76, 82, 124, 128, 130, andl36) and 
those in which they occupy different major shells {A = 
96, 100, 110, 116, 148, 150, 154, 160, andl98). For example 
in 32Ge44 — >34Se42 decay both protons and neutrons oc- 
cupy the shell 28-50, while in 4|"Pd64 — >i|°Cd62 decay 
protons occupy the shell 28-50, and neutrons occupy the 
shell 50-82. The magnitude of the Fermi matrix element 
which is related to the overlap of the proton and neutron 
wave functions is therefore different in these two classes 
of nuclei, being large in the former and small in the latter 
case. This implies a considerable amount of isospin viola- 
tion for nuclei in the first class. The two classes are sepa- 
rated by lines in Table UlT] in order to make the distinction 
clear. For completeness, we have added at the bottom of 
the table the IBM-2 calculation of ''®Ca^'*^Ti, assuming 
""^^Ca to be double magic. IBM-2 is rather poor in this 
evidenced by the large Fermi matrix element, 
and the values in the table for ^^Ca^^^Ti decay should 
be considered a rough estimate. Table HIT] also shows the 



tensor matrix elements, M^?"'. These are systematically 
small (about 5% of M^^'') and have sign opposite to or 

the same as Mq^"^ when protons and neutrons occupy the 
same major shell or not, respectively. This behavior can 
be traced to the fact that the neutrino potential V{r) is 
different for the tensor contribution than for Fermi and 
Gamow- Teller contributions. In the notation of Table 
VIII of Ref. jUl, V{r) = H{r) for Fermi and Gamow- 
Teller matrix elements and V{r) = —rH'ir) for tensor 
matrix elements. 

A point of great interest is the comparison among var- 
ious model calculations of the NME. Up to 2009, the 
methods used were the QRPA and the ISM. In addition 
to these, there are now our calculation (IBM-2) and cal- 
culations based on the density functional theory (DFT). 
Among the QRPA calculations there are those of the 
Tubingen group and those of the Jyvaskyla group. These 
calculations often use different parametrizations of the 
SRC, may or may not include g\ in M(°'^\ use different 
values of qa, and are done in the closure or non-closure 
approximation. Therefore, the comparison among ma- 
trix elements in different models should be taken only 
as relative to a given matrix element, for example, ^^Ge 

In Table HVl we compare our results with those of a par- 
ticular QRPA calculation, QRPA-Tii [Hi and of an ISM 
calculation |1J] with Miller-Spencer (M-S) parametriza- 
tion of the SRC. The IBM-2 and QRPA-Tii results show 
a similar variation with A, while the ISM results are, 
apart from the small value in ^^Ca, clustered around 
~ 2.00 in the entire range A = 76 — 136, and are a fac- 
tor of approximately 2 smaller than results from IBM- 
2 and QRPA-Tii. It should be noted that, due to the 
different approximation made in each model, a range of 
values would be more appropriate. For example, if we 
set to zero the Fermi matrix element in our calculated 
^^Ca ^^^Ti decay, we obtain M^^"'^ = 1.33 and thus our 
matrix elements should be more appropriately quoted as 
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Table III. IBM-2 nuclear matrix elements M'""-* (dimensionless) for decay with Jastrow M-S SRC and gv /g.A ~ 1/1.269. 



A 






14°") 


M(0'')[0+] 








M(0'^)[0+] 


76 

oz 


4.10 

O.ZD 


-2.53 


-0.25 


5.42 


1.81 

U.oD 


-1.21 
-u.Dy 


-0.10 
-u.uo 


2.46 

i.Zo 


96 
100 
110 
116 


2.26 
3.32 
3.22 
2.49 


-0.24 
-0.33 
-0.26 
-0.23 


0.13 
0.20 
0.24 
0.15 


2.53 
3.73 
3.62 
2.78 


0.04 
0.88 
0.41 
0.78 


-0.00 
-0.09 
-0.04 
-0.06 


0.00 
0.05 
0.03 
0.04 


0.04 
0.99 
0.46 
0.85 


124 

128 
130 
136 


Z.D9 

3.46 
3.12 
2.59 


-l.oo 
-1.90 
-1.69 
-1.37 


A 1 O 

-0.13 
-0.16 
-0.14 
-0.11 


3.o0 
4.48 
4.03 
3.33 


O AO 

2.03 
2.44 
2.33 
1.40 


1 oo 

-1.40 
-1.32 
-0.75 


A 1 A 

-0.10 
-0.10 
-0.09 
-0.04 


O '7A 
2. lU 

3.22 
3.07 
1.82 


148 
150 
154 
160 
198 


1.73 
2.03 
2.23 
3.25 
1.64 


-0.28 
-0.28 
-0.26 
-0.31 
-0.23 


0.08 
0.11 
0.12 
0.18 
0.10 


1.98 
2.32 
2.50 
3.62 
1.88 


0.22 
0.35 
0.01 
0.66 
0.07 


-0.04 
-0.05 
-0.01 
-0.08 
-0.01 


0.01 
0.02 
0.01 
0.05 
0.01 


0.25 
0.39 
0.02 
0.75 
0.08 


48 


1.53 


-1.03 


-0.19 


1.98 


3.62 


-3.78 


-0.13 


5.83 



Table IV. Nuclear matrix elements for Qv^fi decay to the 
ground state, 0+, in IBM-2 with Jastrow M-S SRC and 
OA = 1.269, QRPA-Tii with Jastrow M-S SRC and qa = 1.254 
jli, and the ISM with Jastrow M-S SRC and qa = 1.25 [H. 
All matrix elements are in dimensionless units. 



Decay M^°"^ 





IBM-2 


QRPA-Tii 


ISM 




1.98 




0.54 




5.42 


4.68 


2.22 




4.37 


4.17 


2.11 


962i.^96Mo 


2.53 


1.34 




ioomo^iooRu 


3.73 


3.53 




110pd_^110cd 


3.62 






iiecd^iiegn 


2.78 


2.93 






3.50 




2.02 


i28Te_^i28xe 


4.48 


3.77 


2.26 


laorpg^isoxe 


4.03 


3.38 


2.04 


"^Xe^i^^Ba 


3.33 


2.22 


1.70 


i48Nd^i«Sm 


1.98 








2.32 






i^''Sm^i54Q^j 


2.50 






i60Gd^«0Dy 


3.62 






198pt^l98jjg 


1.88 







M(°'') = 1.33-1.98. The sensitivity of the IBM-2 results 
to parameter changes, model assumptions and operator 
assumptions is further discussed in Sec. Ill C 31 The re- 
sults in Table IIVI are summarized in Fig. [S] where they 
are plotted as a function of neutron number. The reason 
for this way of plotting is due to shell effects, as dis- 
cussed in Sec. VI B of Ref. jl5|. The matrix elements 
M^^'') attain their smallest values at the closed proton 
and neutron shells due to the form of the transition op- 
erator which for /3~/3~ decay annihilates a neutron pair 



and creates a proton pair. These shell effects are very 
clear in both the IBM-2 and the QRPA calculations, and 
to some extent also in the ISM calculation. They are re- 
sponsible for the small matrix element in the decay of the 
doubly magic nucleus *^Ca. They are also responsible for 
the ratio of the matrix elements of two different isotopes 
of the same element. For example, a simple calculation 
using the pair operators of Eq. (42) of Ref. [Tsj, gives 
M(0'')(i28Te)/Af(0'''("0Te) = 1.11, to be compared with 
1.11 from IBM-2, 1.13 from QRPA-Tii, and 1.11 from 
ISM. 



Our results to are shown in Table |Vl Because of 
the reduced phase-space factor for decay to Oj, this ta- 
ble is of less interest. In this case, there appears to be 
no correlation between IBM-2 and other calculations. It 
should be noted, however, that the QRPA-Tii results 
shown in Table |V] were done before an error was discov- 
ered in the treatment of short-range correlations (3^ and 
with two different methods for treating the excited 
state, the recoupling method (RCM) and the boson ex- 
pansion method (BEM) jl^. These results are therefore 
inconsistent by a factor of approximately 2 with those 
in Table HV] and Ref. jl^ for 0+ . Also, IBM-2 calcula- 
tions have been done without including intruder configu- 
rations. It is known that, in some cases, the state is an 
intruder state. The most notable cases are Ge, Mo, Cd, 
Nd, and Hg isotopes j3^ p. 180]. Although configuration- 
mixing IBM-2 calculations for these nuclei are available, 
they have not been implemented yet in the calculation 
of to states. The comparison among IBM-2, 

QRPA-Tii, and the ISM is shown in Fig. E] 

The most detailed comparison among different model 
calculations yet has been done recently by Suhonen (36j . 
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Figure 5. (Color online) IBM-2 results for 0!//3/3 compared with QRPA-Tii [Tl| and the ISM [Ti|. 



Table V. Neutrinoless nuclear matrix elements to the first 
excited state, Oj, in IBM-2, QRPA-Tii RCM/BEM and 
the ISM [331. All matrix elements are in dimensionless units. 



Decay 





IBM-2 


QRPA-Tii 


ISM 




5.83 




0.68 


^«Ge^^«Se 


2.46 


1.28/0.99 


1.49 




1.23* 


1.34/0.95* 


0.28* 




0.04 






ioomo^iooRu 


0.99 


1.27/1.76 




iiOpd^iiOCd 


0.46 




ii6Cd^"«Sn 


0.85 






124Sn^l24.pg 


2.70 




0.80 


128Tg_^128Xe 


3.22* 






i30rpg^i30xe 


3.07 




0.19 


"6Xe^i3«Ba 


1.82 


4.42/0.44 


0.49 




0.25 






0.39 






i5*Sm-^i54Q^ 


0.02 






160Qd^l60j3y 


0.75 






198p^_^198jjg 


0.08* 








* Negative C 


J-value 





This author has shown a very close correspondence be- 
tween the IBM-2 results and the QRPA-Jy result. Ta- 
ble I VII and has argued that the reason why the QRPA 
and IBM-2 agree is because the QRPA can be seen as 
a leading-order boson expansion. This statement should 
be taken, however, with caution since the QRPA results 
require the adjustment of the parameter gpp. 
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Figure 6. (Color online) IBM-2 results for Oi/PP decay to 0^ 
compared with QRPA-Tii [H and the ISM [ssl- 



Another question which has been extensively ana- 
lyzed in recent months is the size of the Fermi matrix 
elements, and its comparison among different models. 
To this end, it is convenient to introduce the quantity 
Xf — [gv / gA)'^ Mp'^'^ / Mq^K This quantity is shown in 
Table rvTIl One can see that the situation is more complex 
than in the case of the overall matrix elements, due to the 
different approximations made by the different models. 
The IBM-2 results are large for nuclei in which protons 
and neutrons occupy the same major shells, but small 
in cases in which protons and neutrons occupy differ- 
ent major shells. For QRPA-Jy they are uniformly large 
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Table VI. Comparison between IBM-2 and QRPA-Jy ^ nu- 
clear matrix elements A/''^'^' (dimensionless) for Ou/Sp decay. 



A 0+ Ot 





IBM-2 


QRPA 


IBM-2 


QRPA 


48 


1.98 


1.09-1.89 


5.83 




76 


5.42 


2.28-4.17 


2.46 


2.47-5.38 


82 


4.37 


2.11-3.51 


1.23 


0.831-1.85 


96 


2.53 


2.00-2.07 


0.04 


1.96 


100 


3.73 


2.26-2.74 


0.99 


0.31 


110 


3.62 


3.63-4.51 


0.46 


0.96-1.73 


116 


2.78 


2.36-3.98 


0.85 


0.25 


124 


3.50 


2.58-4.18 


2.70 


3.96-5.88 


128 


4.48 


2.74-4.15 


3.22 




130 


4.03 


2.60-3.78 


3.07 


3.88-6.61 


136 


3.33 


1.83-2.53 


1.82 


2.75-6.08 



(X-F ~ —0.30) while for the ISM they are uniformly small 
iXF - -0.15). 

The large Fermi matrix elements in IBM-2 with pro- 
tons and neutrons in the same major shell and in the 
QRPA throughout point to large isospin violation in the 
wave functions of the initial and final nuclei. In the case 
of 2vj3j3 decay, if isospin is a good quantum number, the 
Fermi matrix elements should identically vanish. By a 
similar argument, the Fermi matrix elements in 0^/3/3 are 
expected to be small, although not zero, the main differ- 
ence between 2vf3/3 and QvjSp being the neutrino poten- 
tial, given in Appendix A, Table IXIXI For this reason, 
the calculated values of xf given in Table IVIII may be 
entirely spurious and should be considered with a large 
error. It is difficult to estimate the error in xf introduced 
by isospin violation, since there are no direct experimen- 
tal data for single-/? 0+ — s> 0+ transitions from odd-odd 
to even-even (or vice versa) nuclei in heavy nuclei. The 
estimate of the error also depends on which nucleus is 
considered. The maximum error is 100%, if the Fermi 
matrix elements are entirely spurious, in which case the 
values of Table |VIl] should be quoted as -0.42(42) for 
''^Ca and similarly for all others. Another estimate is to 
extract the error by comparing with the ISM calculations 
which have the smallest values of x^?, in which case the 
quoted value should be —0.42(28) with an error of 67%. 
We have used this estimate of the error in the following 
subsection 3 and in Table Kill 



Table VII. Comparison among Fermi matrix elements, xf, in 
IBM-2, QRPA-Jy [H and the ISM [Hi. 



Decay 




Xf 




IBM-2 


QRPA-Jy 


ISM 


■^^Ca 


-0.42 


-0.56^ 


-0.14 


'^''Ge 


-0.38 


-O.22" 


-0.10 


**^Se 


-0.42 


-0.28 


-0.10 




-0.06 


-0.43 




i°°Mo 


-0.06 


-0.40 




llOp^j 


-0.05 


-0.38 


-0.16 


"^Cd 


-0.06 


-0.28 


-0.19 


'^''Sn 


-0.35 


-0.42 


-0.13 


i28Te 


-0.34 


-0.37 


-0.13 


""Te 


-0.34 


-0.37 


-0.13 


"6Xe 


-0.33 


-0.34 


-0.13 


"«Nd 


-0.10 






150Nd 


-0.09 






i5*Sm 


-0.07 






i«°Gd 


-0.06 






198p^ 


-0.09 







Reference [H. 



2. QvPP decay with heavy neutrino exchange 



These matrix elements can be simply calculated by re- 
placing the potential v{p) = 2t:~'^[p\p + A)\~'^ of Eq. ([13]) 
with the potential Vh{p) = 27r^^(memp)^^ of Eq. (|14p . 
Table IVIlTl gives the corresponding matrix elements. The 
index h distinguishes these matrix elements from those 
with light neutrino exchange. Our results are compared 
with QRPA-Tii results in Table HXl These QRPA results 
are obtained with Jastrow SRC [3j| and prior to the more 
refined treatment of SRC of Ref . [131 , and are shown here 
for the sake of comparison of the A dependence, not of 
their absolute magnitude, which is a factor of approxi- 
mately 2 smaller than in IBM-2. It has been suggested 
that measurement of neutrinoless double-/? decay in dif- 
ferent nuclei may be used to distinguish between the two 
mechanisms, light or heavy neutrino exchange. However, 
the results in Table IIIII and Table IIXI are highly corre- 
lated as is clear from the fact that they are obtained one 
from the other by replacing the potential v{p) with Vh[p). 
Therefore, this criterion cannot be used to distinguish be- 
tween the two mechanisms 1451 . 



In addition to the calculations discussed above, sev- 
eral others have been made, most notably in the de- 
formed QRPA (40l,|4l| and in the projected Hartree-Fock- 
Bogoliubov (HFB) framework [43 , and using the energy 
density functional method jlsj . Since these use SRC with 
Argonne/CD-Bonn and the unitary correlation method 
(UCOM) they will be discussed at the end of Sec. HTCSl 



3. Sensitivity to parameter changes, model assumptions 
and operator assumptions 

Many ingredients go into the calculation of the nuclear 
matrix elements. In Ref. fisj . the sensitivity to input 
parameter changes was estimated from the sensitivity to 
parameter changes in five quantities: (1) single-particle 
energies; (2) strengths of interactions: (3) oscillator pa- 
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Table VIII. Nuclear matrix elements for the heavy neutrino exchange mode of the neutrinoless double-/? decay to the ground 
state (columns 2, 3, 4, and 5) and to the first excited state (columns 6, 7, 8, and 9) using the microscopic interacting boson 
model (IBM-2) with Jastrow M-S SRC and gv/gA = 1/1.269. 



A 0+ 0^ 





Oil, 


I'h 


J J, 


A/,r^[0+] 






J-h 


a4°-)[o+1 


76 
82 


52.5 
43.9 


-39.5 
-33.9 


-29.0 
-29.4 


48.1 
35.6 


19.5 
8.53 


-15.4 
-7.36 


-10.5 
-4.96 


18.5 
8.14 


96 
100 
110 
1 1 fi 


33.6 
54.5 
49.3 


-18.8 
-30.0 
-26.5 
9n 1 


13.6 
26.1 
29.9 


59.0 
99.3 
95.7 

fi7 1 


0.700 
4.59 
18.4 

1 Q 


-0.382 
-2.55 
-9.73 

7 79 


0.321 
1.70 
10.0 

7 07 


1.26 
7.87 
34.5 

9^ 8 


124 
128 
130 
136 


37.0 
47.0 
42.7 
33.7 


-27.1 
-34.1 
-30.9 
-24.3 


-16.1 
-19.7 
-17.9 
-13.7 


37.8 
48.4 
44.0 
35.1 


25.2 
28.0 
25.9 
12.9 


-19.0 
-21.1 
-19.6 
-9.90 


-10.6 
-11.1 
-10.0 
-4.39 


26.4 
30.0 
28.1 
14.7 


148 
150 
154 
160 
198 


35.8 
39.9 
38.6 
51.9 
34.7 


-21.2 
-23.0 
-21.8 
-28.7 
-20.8 


10.4 
14.2 
14.9 
23.2 
13.8 


59.4 
68.4 
67.1 
92.9 
61.5 


3.99 
5.90 
1.93 
13.1 
1.51 


-2.43 
-3.42 
-1.12 
-7.34 
-0.899 


1.07 
1.90 
1.20 
6.65 
0.688 


6.57 
9.92 
3.82 
24.3 
2.76 


48 


26.8 


-20.1 


-22.9 


16.3 


23.8 


-27.3 


-6.95 


33.8 



Table IX. Neutrinoless nuclear matrix elements to the ground 
state, Ot, in IBM-2 and QRPA-Tii [1| for the heavy neutrino 
exchange mode and M-S SRC. All matrix elements are in 
dimensionless units. 



Decay 




Mr 


IBM-2 


QRPA-Tii 


48ca^«Ti 


16.3 






48.1 


32.6 


s^Se^^^Kr 


35.6 


30.0 


seZr^'^^Mo 


59.0 


14.7 


iooMo^iooRu 


99.3 


29.7 




95.7 




ii6cd^"6Sn 


67.1 


21.5 




37.8 




i28Te^i28xe 


48.4 


26.6 


lao^g^iaoxe 


44.0 


23.1 


"^Xe^i^'^Ba 


35.1 


14.1 


i4»Nd^"«Sm 


59.4 




«°Nd^i«'Sm 


68.4 


35.6 


i^^Sm^is^Gd 


67.1 




160gj^l60j)y 


92.9 




198pt^l98jjg 


61.5 





rameter in the single-particle wave functions; (4) closure 
energy in the neutrino potential: and (5) nuclear radius. 
(1) The sensitivity to single-particle energies was empha- 
sized in Ref. 1^461 within the framework of the QRPA 
and has been the subject of several experimental inves- 
tigations [13, We estimate it to be 10%. (2) We 
estimate the sensitivity strengths of interactions, in the 
present case the surface delta interaction used to calcu- 
late the structure of the pair states. We estimate this to 



be 5%. We note that this is the main source of sensitivity 
in the QRPA, especially in nuclei close to the spherical- 
deformed transition, for example ^^°Nd. (3) We estimate 
the sensitivity to the oscillator parameer to be 5%. (4) 
We estimate the sensitivity to closure energy to be 5%. 
(5) If the matrix elements are quoted in dimensionless 
units there is also the sensitivity to R. We estimate this 
to be 5%. However, this sensitivity can be reduced to 
a small value, 1%, if the experimental rms value is used 
instead of the formula R = RqA^^^. The total estimated 
sensitivity to input parameters is 30% if all contributions 
are added or 14% if combined in quadrature. 

In addition, we estimate the sensitivity to model as- 
sumptions to be: (1) truncation to pairs with angular 
momentum J = and J = 2 (S-D space) and (2) isospin 
purity. We estimate the former to range from 1% in 
spherical nuclei to 10% in deformed nuclei. For the lat- 
ter we estimate this to be small, 1%, for GT and T ma- 
trix elements, and large, 10%, for F matrix elements. 
Taking into account the fact that the F matrix elements 
contribute only ~20% to the total matrix elements, we 
estimate the total sensitivity to model assumptions to 
be ranging from 3% in spherical nuclei, to 12% in de- 
formed nuclei (in addition) or 2%-10% (in quadrature). 
A special case is that of '^^Ca decay for which, for reasons 
mentioned above, the Fermi matrix elements are overesti- 
mated. In this case the sensitivity to model assumptions 
may be as high as 20% (in addition) or 16% (in quadra- 
ture) . 

Finally, there is the estimated sensitivity to operator 
assumptions: (1) Form of the transition operator. We 
have already commented in Sec. II. A on the differences 
between Simkovic's Q and Tomoda's Q formulations. 
This is a source of considerable uncertainty. We esti- 
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Table X. Comparison among matrix elements M'-"''' calculated with Miller-Spencer (M-S) and Argonne/ CD-Bonn (CCM) in 
IBM-2, with M-S, CCM, and UCOM in QRPA-Tii [H, [H, HJl , and with M-S and UCOM in the ISM [li,[3l|. Note that the 
QRPA matrix elements are evaluated using qa ~ 1.254 and the ISM matrix elements are evaluated using qa = 1.25. 



A IBM-2 QRPA-Tii ISM 





M-S 


CCM 


M-S'' 


ccm'^ 


TTCOM'' 


M-S"^ 


TICnM'* 


48 




1.98 


2.28/2.38 








0.59 


0.85 


76 




5.42 


5.98/6.16 


4.68 


5.81/6.32 


5.73 


2.22 


2.81 


82 




4.37 


4.84/4.99 


4.17 


5.19/5.65 


5.09 


2.11 


2.64 


96 




2.53 


2.89/3.00 


1.34 


1.90/2.09 


1.79 






100 




3.73 


4.31/4.50 


3.52 


4.75/5.25 


4.58 






1 1 n 




3.62 


'^.10/ 41,01 












116 




2.78 


3.16/3.29 


2.93 


3.54/3.99 


3.54 






124 




3.50 


3.89/4.02 








2.02 


2.62 


128 




4.48 


4.97/5.13 


3.77 


4.93/5.49 


4.76 


2.26 


2.88 


130 




4.03 


4.47/4.61 


3.38 


4.37/4.92 


4.26 


2.04 


2.65 


136 




3.33 


3.67/3.79 


2.22 


2.78/3.11 


2.76 


1.70 


2.19 


148 




1.98 


2.36/2.49 












150 




2.32 


2.74/2.88 




3.34^ 








154 




2.50 


2.91/3.04 












160 




3.62 


4.17/4.34 




3.76^ 








198 




1.88 


2.25/2.37 












" Reference 
^ Reference 
Reference 
Reference 
Reference 


2e 

14 
33 
41 


I- 
I- 



mate the sensitivity to 5% by comparing our calculations 
using Tomoda's and Simkovic's formulations. However, 
there still remains the question of the recoil contribution 
which was a source of major disagreement in early calcu- 
lations. (2) Finite nuclear size (FNS). We estimate this 
to be small, 1%, due to the fact that we have used re- 
alistic nucleon form factors with parameters determined 
from experiment. (3) Short-range correlations (SRC). 
The sensitivity to the form of short-range correlations has 
been very recently the subject of many studies. To in- 
vestigate this point, we have done calculations with three 
different types of correlations, (a) Jastrow Miller-Spencer 
(MS) and (b/c) Argonne/CD-Bonn (CCM). The results 
are shown in Table [X] where we also show a comparison 
with results of calculations in the QRPA and the ISM 
using M-S, CCM and the unitary correlation operator 
method (UCOM) jH]. It appears that in going from M- 
S to CCM or UCOM the matrix elements in all three 
methods (IBM-2, the QRPA, and the ISM) increase. In 
IBM-2 the multiplicative factor when going from M-S to 
CCM-Argonne is 1.10-1.20. In QRPA-Tii it is 1.21-1.42 
from M-S to CCM-Argonne and 1.21-1.33 from M-S to 
UCOM. In the ISM the factor is 1.25-1.30 from M-S to 
UCOM. This multiplicative factor was taken into account 
in Ref. j50| when comparing IBM-2 matrix elements with 
the ISM. The discrepancy between IBM-2 and QRPA-Tii 
multiplicative factors is not understood and should be in- 
vestigated further. In Table 1X1 we have also added recent 
calculations in the deformed QRPA for decay of ^^"Nd 
and ^^''Gd [4l|. One may note that the correspondence 



between the QRPA and IBM-2 persists even in deformed 
nuclei. 

The total estimated sensitivity here is 11% (addition) 
or 7% (quadrature), under the assumption of no recoil 
contribution to the matrix elements. Combining all con- 
tributions, we have a total estimated sensitivity of 44%- 
55% if all the contributions are added or 16%-19% if they 
are combined in quadrature. 

The short-range correlations affect decay differ- 

ently than Oiy/3/3. We have therefore investigated the de- 
pendence of QvhPP matrix elements with M-S, and CCM 
correlations. The results are shown in Table IXII We see 
here an increase of a factor from 1.69 to 2.80 when going 
from M-S to Argonne/CD-Bonn. This is because, as re- 
marked in Appendix A, the neutrino potential for heavy 
neutrino exchange is a contact interaction in configura- 
tion space and thus strongly influenced by SRC. The cor- 
relation function in Eq. (|17p has a value /./(O) = 1 — c 
at r = 0. For the M-S parametrization c = 1, /./(G) = 0, 
and thus, in the absence of a nucleon form factor, the 
matrix elements M^^^ vanish. The Argonne/CD-Bonn 
parametrizations have c = 0.92/0.46 and thus a non-zero 
value at r = 0. These results are modified by the pres- 
ence of the nucleon form factors of Eq. ((HI), and the 
final results depend strongly on the choice of gv (j>^ ) and 
From Table IXll columns 4 and 5, it appears that 
the increase in the matrix elements when going from M-S 
to CCM-Argonne in QRPA-Tii is much larger, from 7.01 
to 10.1, than in IBM-2. As in the case of light neutrino 
exchange, this discrepancy is not understood and should 
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be investigated further. The large increase both in IBM-2 
and the QRPA also points to the strong sensitivity of the 
calculated QvhPP matrix elements to the specific form of 
SRC, and thus to the fact that the treatment here and in 
other calculations in the literature, through the nucleon 
form factors, may not be satisfactory. A more consis- 
tent treatment is discussed in Refs. fsil Is^. In view of 
all these problems, we estimate the sensitivity to SRC for 
0^/?l/3/3 decay to be much larger (50%) than that for Qvj3(3 
decay (5%). The total estimated sensitivity to operator 
assumptions for QvhPfi is then 56% (addition) and 50% 
(quadrature). 



Table XL Comparison among m'^'^^ matrix elements calcu- 
lated with different Jastrow parametrizations for the SRC: 
Miller-Spencer and CCM (Argonne and CD-Bonn) in IBM-2 
and QRPA-Tii 



A 




IBM-2 




QRPA-Tii 


M-S 


CCM 


M-S'' 


CCM'' 


48 


16.3 


46.3/76.0 






76 


48.1 


107/163 


32.6 


233/351 


82 


35.6 


84.4/132 


30.0 


226/340 


96 


59.0 


99.0/135 


14.7 




100 


99.3 


165/224 


29.7 


250/388 


110 


95.7 


155/208 






116 


67.1 


110/149 


21.5 




124 


37.8 


79.6/120 






128 


48.4 


101/152 


26.6 




130 


44.0 


92.0/138 


23.1 


234/364 


136 


35.1 


72.8/109 


14.1 




148 


59.4 


103/142 






150 


68.4 


116/160 


35.6 




154 


67.1 


113/155 






160 


92.9 


155/211 






198 


61.5 


104/141 







^ Reference 
^ Reference 



To summarize the situation we show in Table IXIII our 
final results with M-S SRC, together with an estimate of 
the error, based on the arguments given above. The error 
estimate is 30% in ^^Ca, 19% in nuclei with protons and 
neutrons in the same major shell and 16% in nuclei with 
protons and neutrons in different major shells, for Oz^/3/3. 
For Oi/fi/Sp, our estimated error is dominated by SRC. In 
Table En] we have used 58% in ''^^Ca, 53% in nuclei with 
protons and neutrons in the same major shell and 52% 
in nuclei with protons and neutrons in different major 
shells. 

Finally, having investigated the effect of short-range 
correlations on 0iy/3(3 we are now able to compare our 
results with all available calculations done with the same 
SRC including DFT jHI and HFB 113. These are shown 
in Fig. [3 We note now that while the ISM/QRPA/IBM- 
2 have the same trend with A, the other two do not. For 
the isotopic ratio M'^^'^^^^^Te)/ M^°''^^^°Te) the DFT 
method gives 0.86 in sharp contrast with the value 1.11. 



Table XII. Final IBM-2 matrix elements with M-S SRC and 
error estimate. 

Decay Light neutrino exchange Heavy neutrino exchange 



*«Ca 
^^Ge 
^^Se 
9«Zr 
Mo 
"Opd 
"«Cd 
i^*Sn 
i28Te 
"OTe 
i36Xe 
"«Nd 
"ONd 
i^*Sm 
i«OGd 
««Pt 



1.98(59) 
5.42(103) 
4.37(83) 
2.53(40) 
3.73(60) 
3.62(58) 
2.78(44) 
3.50(67) 
4.48(85) 
4.03(77) 
3.33(63) 
1.98(32) 
2.32(37) 
2.50(40) 
3.62(58) 
1.88(30) 



16.3(95) 
48.1(255) 
35.6(189) 
59.0(309) 
99.3(516) 
95.7(498) 
67.1(321) 
37.8(200) 
48.4(257) 
44.0(233) 
35.1(186) 
59.4(309) 
68.4(356) 
67.1(349) 
92.9(483) 
61.5(320) 



Also, while the ISM/QRPA/IBM-2 have a small value 
for ^^Zr, DFT has a large value. We therefore conclude 
that the approximations made in the DFT/HFB lead to 
a different behavior with A. This point is currently being 
investigated jlj. Also, the Fermi matrix elements in the 
DFT are comparable to those in IBM-2 and larger than 
those in the ISM El. 




i • IBM-2 

1 ▲ QRPA-Tii 

*♦ ^ ^ i T QRPA-Jy 

* Qd r : EDF 

T ♦♦*^Nj,f Pt 
i * t * 



60 80 
Neutron number 



ICQ 



120 



Figure 7. (Color online) IBM-2 results for OvPfi nuclear ma- 
trix elements compared with QRPA-Tii [11 . the ISM [l3 |. 
QRPA^ [H, [HiSi, QRPA-deformed [JJ DFT and 
HFB [ig. 
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D. Limits on neutrino mass 

1. Light neutrino exchange 

The calculation of nuclear matrix elements in IBM-2 
can now be combined with the phase-space factors cal- 
culated in and given in Table III and Fig. 8 of that 
reference to produce our final results for half-lives for 
light neutrino exchange in Table IXIIII and Fig. HI The 
half-lives are calculated using the formula 



(21) 



We note here that the combination must be done con- 
sistently. If the value of gA is included in Mqi^, then it 
should not be included in g'^qJ and similarly for a factor of 

4 included in some definition of G^^J and not in others 
• See Eq. (53) of Ref. This point has caused con- 
siderable confusion in the literature. In Table IXIIII and 
Fig. [5] the values (mi,) = 1 eV and gA = 1.269 are used. 
For other values they can be scaled with \ {m^) /me| and 
9\- 



Table XIII. Left: Calculated half-lives in IBM-2 M-S SRC for 
neutrinoless double- /3 decay for (m„) = 1 eV and gA = 1.269. 
Right: Upper limit on neutrino mass from current experi- 
mental limit from a compilation of Barabash [s^ . The values 
reported by Klapdor-Kleingrothaus et al. |59l |. the IGEX Col- 
laboration (60I I , a nd the recent limits from KamLAND-Zen 
[61II and EXO mal are also included. 



Decay 


Of /I n24 N 

^1/2(10 yr) 




(m„> (eV) 




1.03 


> 5.8 X 10^2 


< 4.2 


'"^Ge^^'^Se 


1.45 


> 1.9 X 10^^ 


< 0.28 






1.2 X 10^^^ 


0.35 






> 1.6 X 10~^ 


< 0.30 


*^Se^*^Kr 


0.52 


> 3.6 X 10^~ 


< 1.2 




0.77 


> 9.2 X 10^^ 


< 9.1 


100mo->i«0Ru 


0.46 


> 1.1 X 10^'' 


< 0.64 


iiOpd_,"Ocd 


1.60 






ii6Cd^ii6Sn 


0.78 


> 1.7 X 10^^ 


< 2.1 




0.91 






i28Te^i28xe 


8.53 


> 1.5 X lO^"* 


< 2.4 


I30^g^i30xe 


0.44 


> 2.8 X 10^'' 


< 0.39 


i36xe^"«Ba 


0.62 


> 5.7 X 10^^^ 


< 0.33 






> 1.6 X 10^^3 


< 0.20 




2.54 








0.30 


> 1.8 X 10^^ 


< 4.1 




5.34 








0.80 






198p^_^198jjg 


3.77 







Reference 
Reference 
Reference 
Reference 

The effective neutrino mass is the quantity we want 
to extract from experiment. Unfortunately, the axial 




0.2L 



50 



Mo 

Nd 

100 150 
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200 



Figure 8. (Color online) Expected half-lives for (m^) = 1 eV, 
gA = 1.269. The points for ^^^Te and "**Nd decays are not 
included in this figure. The figure is in semilogarithmic scale. 



vector coupling constant is renormalized in nuclei to 
gA,eff- A (model-dependent) estimate of gA,eff can be 
obtained from the experimental knowledge of single-/? de- 
cay and/or of 2i//3/3 decay. This will be discussed in the 
following section. Here we show in Fig. |9] and Table IXlII] 
the limits on neutrino mass from current experimental 
upper limits using IBM-2 matrix elements of Table IVl 
and gA = 1.269. In addition to the experimental upper 
limits, a value has been reported for the half-life in ^^Ge, 
1.2 X lO^Vr jUl. This is also reported in Fig. El 

The average light neutrino mass is constrained by at- 
mospheric, solar, reactor and accelerator neutrino oscil- 
lation experiments to be (Tsj 



{nil,) = \cl2cl2m1 -1-0^35^2^26"^^ 



cos 



/ 2 2 2\ '"l 

(mi, 7712,7713) = 



sin-dij, .^2,3 = [0, 27r] 



(22) 



Srn^ 5m? 



, ±Am' 



Using the best fit values (75| 



sin^ t?i2 = 0.213, sin^ i?i3 



0.016, 



sin^ i?23 = 0.466, dm^ = 7.67 x 10"^ eV^ 
Am^ = 2.39 X 10"^ eV^ 



(23) 



we obtain the values given in Fig. 1101 In this fig- 
ure we have added the current limits, for a a — 1.269, 
coming from CUORICINO [H, IGEX Eg, NEMO-3 
H, KamLAND-Zen [6l|, and the EXO jHll experiment. 
Also, henceworth we use c = 1 to conform with standard 
notation. 
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Figure 9. (Color online) Limits on neutrino mass from current 
experimental limits from a compilation of A. Barabash (58| 
and recent measurement for ^^^Xe from EXO [6^ . The values 
reported by Klapdor-Kleingrothaus et al. (59| is shown by the 
symbol X. The figure is in semilogarithmic scale. The shaded 
area represents the values of allowed by the current 

experiments. 
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Figure 10. (Color online) Current limits to (rriv) from 
CUORICINO IpL IGEX [60|, NEMO-3 [H, KamLAND-Zen 
[H, and EXO [iil, and IBM-2 M-S SRC nuclear matrix ele- 
ments. The value of Ref. [59l | is shown by an X. It is consistent 
only with nearly degenerate neutrino masses. The figure is in 
logarithmic scale. 



2. Heavy neutrino exchange 

The half-hves for this case are calculated using the for- 
mula 



'^1/2 1 



= G("j|Afo.j'|7?f 



7] = 



(24) 



k— heavy 



The expected half- lives for \ri\ = 2.75 x 10"^, and using 
the IBM-2 matrix elements of Table HXl are shown in Ta- 
ble |XIVl For other values of -q they scale as |r/p. There 
are no direct experimental bounds on rj. Recently, Tello 
et al. [7j have argued that from lepton flavor violating 
processes and from Large Hadron Collider (LHC) exper- 
iments one can put some bounds on the right-handed 
leptonic mixing matrix Uek, heavy and thus on rj. In the 
model of Ref. (7(, when converted to our notation, i] can 
be written as 



E 

-heavy 



.2 rrip 



(25) 



where M\y is the mass of the M^-boson, Mw = (80.41 ± 
0.10) GeV [2I, MwR is the mass of the WR-hoson, 
assumed in Q to be Mwr = 3.5 TeV and V = 
{MwR/Mwf U. The ratio {Mw/MwRf is then 2.75 x 



10^^, the value we have used in the left portion of Ta- 
ble IXIVI By comparing the calculated half-lives with 
their current experimental limits, we can set limits on the 
lepton-nonconserving parameter I77I, shown in Table IXTVl 
and Fig. [Til 



If we write 



V = 



M^R 



(26) 



we can also set limits on the average heavy neutrino mass, 
{niu^ ) , as shown in the last column of Table IXIVI This 
limit is model dependent since it is assumed that Mwr = 
3.5 TeV. For other values of Mwr it scales as M^^. 

If both light and heavy neutrino exchange contribute, 
the half-lives are given by 



'1/2J 



a 



(0) 



(27) 



It is interesting to note here the possibility of interference 
between light and heavy neutrino exchange, as empha- 
sized recently by several authors. The limits presented 
in Figs. [9] [101 and El £^re based on the calculation with 
M-S SRC. If CCM SRC are used, they should be multi- 
plied by ^ 1.2 (light neutrino exchange) and ^ 2.0 (heavy 
neutrino exchange). 
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Table XIV. Left: Calculated half-lives for neutrinoless double-/? decay with exchange of heavy neutrinos for rj = 2.75 x 10^^ and 
QA = 1.269. Right: Upper limits of |Ti| and lower limits of heavy neutrino mass (see text for details) from current experimental 
limit from a compilation of Barabash |5^ . The value reported by Klapdor-Kleingrothaus et al. [59l |. and the IGEX Collaboration 
[60l |. and the recent limits from KamLAND-Zen [gH and EXO [63 | are also included. 



Decay 


^1/2 (.I'J 




1^1 M n— 6^ 


(m^^)(GeV) 


48ca^4«Ti 


0.77 


> 5.8 X 10^^ 


<1.00 


>0.26 




0.95 


> 1.9 X lO^'^ 


<0.061 


>4.2 






1.2 X 10^^^ 


0.077 


3.362 






> 1.6 X 10 


<0.066 


>3.88 


82se^»2j^r 


0.40 


> 3.6 X 10^" 


<0.29 


>0.89 


^^Zr^^'^Mo 


0.07 


> 9.2 X 10^^ 


<0.77 


>0.34 




0.03 


> 1.1 X 10^* 


<0.0047 


>5.5 


iiopd^iiocd 


0.12 








116cd^ll6gn 


0.07 


> 1.7 X 10^^ 


<0.17 


>1.5 


Sn-i^^^'^Te 


0.39 








i28Te^i28xe 


3.71 


> 1.5 X 10^* 


<0.43 


>0.60 


isorpe^iaoxe 


0.19 


> 2.8 X 10^* 


<0.071 


>3.6 


"6Xe^"6gg^ 


0.29 


> 5.7 X lO^-*^ 


<0.061 


>4.2 






> 1.6 X 10^^3 


<0.116 


>2.2 


i48Nd^"8Sm 


0.14 








isoNd^isogni 


0.02 


> 1.8 X 10^^ 


<0.27 


>0.96 




0.38 








leoGd^ieoDy 


0.06 








198p^^l98jjg 


0.18 









^ Reference 5£ | . 

^ Reference 6C | . 

Reference 61 1 . 

Reference 62 1 . 
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Figure 11. (Color online) Limits on the lepton nonconserving 
parameter |77|. The value of Ref. [s^ is shown by an X. The 
figure is in semilogarithmic scale. The shaded area represents 
the values of I77I allowed by the current experiments. 



III. 2uPl3 DECAY 

A. Matrix elements 

As mentioned in the previous section, the calculated 
matrix elements of the GT operator in single-/? decay 
appear to be systematically larger than those extracted 
from the measured ft values of allowed GT transitions. 
To take into account these results, it has been found con- 
venient to renormalize the value of gA to be used in a 
particular model calculation by introducing an effective 
gA,eff.p defined as 



9A,eff,l3 
9A 



\M, 



exp^j3 I 



(28) 



where gA = 1.269 and Mp are the matrix elements for 
single-/? decay. The ratio {gA,eff /9a) is also called the 
quenching, q, or hindrance, h = 1/q, factor. The quench- 
ing of gA comes from two effects: (i) the limited model 
space in which the calculation is done and (ii) the con- 
tribution of non-nucleonic degrees of freedom. A,.... The 
first type of quenching depends strongly on the size of the 
model space used in the calculation and is thus model de- 
pendent. It was extensively investigated in light nuclei, 
A ~ 20, in the 1970s |65l - [67| within the framework of the 
ISM where it was found that gA,eff — 1-0, q = 0.75. In 
heavy nuclei, of particular interest in this paper, the ques- 
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tion of quenching was first discussed by Fujita and Ikeda 
[68| in 1965. These authors analyzed /3-decay in mass 
A ~ 120 nuclei within the framework of various models 
(pairing, pairing plus quadrupole, etc.) and found very 
small quenching factors, q ~ 0.2 — 0.3 thus stimulating 
the statement that massive renormalization of gA occurs 
in heavy nuclei (67) . The second type of quenching was 
extensively investigated theoretically in the 1970s (69l - 
l7l| . This effect does not depend much on the nuclear 
model used in the calculation, but rather on the mecha- 
nism of coupling to non-nucleonic degrees of freedom. It 
is being re-investigated currently within the framework of 
chiral effective field theory (EFT) j72| and there are hints 
that it has a complex structure, in particular that it may 
depend on momentum transfer and that it may lead in 
some cases to an enhancement rather than a quenching. 

The values of gA.eff depend crucially on the model 
used through the size and composition of the model 
space, especially on whether or not spin-orbit partners 
are included in the calculation. For example, while the 
QRPA includes spin-orbit partners, IBM-2 and the ISM 
do not. Conversely, in the ISM calculations the size of 
the model space is 10^, while in the QRPA and IBM- 
2 it is much smaller. In order to extract gA,eff,p for 
a given mass number A, one needs to do a calculation 
of single-/? decay in that region and compare with ex- 
periment where available. Within the context of IBM-2, 
some calculations were done in the 1980s [t^. Very re- 
cently, the problem has been readdressed and results will 
be published soon 

Double- /3 decay depends on gA as g\ and thus its 
quenching is of extreme importance. Since 2i'l3p de- 
cay has now been measured in several nuclei, it pro- 
vides another way to estimate gA,eff which we denote 
by gA.ef f,2ui3p- In this section, we attempt an estimation 

9^A^f~f2vi3p within the framework of IBM-2 in the clo- 
sure approximation and also extract g^A^f f 2vp0 within 
the framework of the ISM in the non-closure approxi- 
mation. The extraction of g^A^ff\vpp non-closure 
approximation will be presented in the forthcoming pub- 
lication mentioned above [t^ . 

2i//3/3 is a process allowed by the standard model 
and thus in principle exactly calculable. The theory of 
2i/p[3 decay was developed by Primakoff and Rosen 
Konopinski (t^, Doi et al. p| and Haxton et al. 
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The calculation of 2vjij3 turns out to be more complex 
than that of 0iyf3/3. 

(i) The closure approximation may not be good and 
one needs to evaluate explicitly the matrix elements to 
and from the individual 1^ and 0^ states in the inter- 
mediate odd-odd nucleus. Fig. [T^ 



0/) 



and 



M 



{2u) 
F,N 



2mec'^) + E,+ - Ei 



0^)R||rt||0^) 



{Qfifi + 2TOeC2) + Eg+ - Ei 



(29) 
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Figure 12. The decay i^^Mosg -^IfRnse, an example of 2^/313 
decay. 



This evaluation has been done in selected nuclei within 
the framework of the pnQRPA [t^, the proton-neutron 
microscopic anharmonic vibrator approach (pnMAVA) 
Izi, [zil, and the ISM jH, H [III and it has been 
programmed very recently within the framework of the 
proton-neutron interacting boson-fermion model (IBFM- 
2) [tII. The calculation requires the difficult task of de- 
termining the structure of the intermediate odd-odd nu- 
cleus. 

(ii) The PSFs cannot be exactly separated from the 
nuclear matrix elements. These factors must be calcu- 
lated separately for each state 1^ and 0^. In order to 
calculate half-lives and other observable quantities the 

product of the PSFs jy ^^^^ matrix elements must be 
calculated and the contributions summed over all indi- 
vidual states. This is a daunting problem compounded 
by the fact that in most calculations the giant Gamow- 
Teller resonance that contributes to the matrix elements 
is not included in the model space. 

The separation of PSFs and nuclear matrix elements 
can be done in two cases: (1) the closure approximation 
(CA) an d (2) t he single-state dominance (SSD) approx- 
imation [82h86| . In both cases, the inverse half-life can 
be written as 



'1/2 



= G 



(0) 



(31) 



In the CA, the matrix elements Moi, can be written as 



M2,^g\M^^-\ 



2v) 



M, 



{2u) 
GT 



A 



GT 



9V_ 

9A 



2 ^(2.) 



Af 



(32) 
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where 



M, 



GT 



E 



M 



{2u) 



0+ 



/ J ' n' n' 



(33) 



The closure energies Act and Ap are defined by 



E, 



Af 



2 {Qpp 



2meC^) + {Ea+^N)-Ei, 



(34) 



where (Ej^) is a suitable chosen excitation energy in the 
intermediate odd-odd nucleus. The Fermi matrix ele- 
ments are suppressed by isospin considerations and are 
often neglected. In case they are not, care must be taken 
since the closure energy Ap is difi^erent from the closure 
energy Act, although, for simphcity Agt ~ Ap = A is 
often used. 

In the SSD approximation, the matrix elements are 
given by 



M, 



{2u) 

GT,SSD 



M 



{2u) 
F,SSD 



(0+||rta||l+)(l+||rta||0+) 
i {Qpp + 2mec2) + E^+ - Ei ' 

(0+||rt||0+)(0+||rti|0+) 
^{Qpp + 2meC^) + E^^ ~Ei' 



(35) 



where E,+ and i?n+ are the energies in the intermediate 
odd-odd nucleus of the single state that dominates the 
decay. From these one can form the quantities 



M, 



GT.SSD 



SSD 



9v \ .A2u) 



and calculate the half- lives from (25), with Gj*^^ given by 

^2v,SSD- 



(36) 



B. Results 

In this article, we present results of a calculation of 
the nuclear matrix elements for 2z//3/3 in the CA using 
the transition operator of Sec. II. A. In this case only 
the terms hyy and /i^J are considered. An advantage 
of the closure approximation for 2j//3^ decay is that the 
nuclear matrix elements can be calculated using the same 
method discussed in Sec. II, by simply replacing the 
neutrino potential v{p) by 



V2v{p) 



5{p) 



(37) 



which is the Fourier-Bessel transform of the configuration 
space potential V{r) = 1. Since our purpose here is a 



Table XV. 2vPj5 matrix elements (dimensionless) to the 
ground state (columns 2 and 3) and to the first excited state 
(columns 4 and 5) using the microscopic interacting boson 
model (IBM-2) in the closure approximation. 



A 




0^ 
















76 
82 


4.34 
3.50 


-2.69 
-2.39 


-1.35 
-0.83 


1.99 
1.03 


96 
100 
110 
116 


2.22 
2.94 
2.98 
2.31 


0.02 
0.03 
0.03 
0.02 


0.00 
0.00 
0.01 
0.01 


0.04 
0.39 
1.48 
0.86 


124 
128 
130 
136 


2.80 
3.63 
3.31 
2.76 


-1.60 
-1.99 
-1.79 
-1.44 


-1.34 
-1.53 
-1.45 
-0.83 


2.15 
2.65 
2.59 
1.63 


148 
150 
154 
160 
198 


1.24 
1.54 
1.91 
2.99 
1.00 


0.02 
0.02 
0.02 
0.02 
0.01 


0.00 
0.00 
-0.00 
0.01 
0.00 


0.17 
0.29 
0.05 
0.51 
0.03 


48 


1.57 


-1.08 


-4.77 


5.02 



direct comparison of ivjSp and QvjBp decays, this avoids 
possible systematic and accidental errors. The CA is not 
expected to be good for 2i^/3/? decay, since only 1+ and 
0"^ intermediate states in the odd-odd nucleus contribute 
to the decay. We use it here only as an estimate, with 
appropriately chosen closure energy 0, p.71] in order to 
extract gA.ef f,2vi3p = gA,eff, which is the purpose of this 
section. 

Our calculated matrix elements for 2i//3/3 decay are 
shown in Table IXVI Also here as in Table IIIII we have 
separated the class of nuclei where protons and neutrons 
occupy the same major shell from those where they do 
not, and we added A = 48 at the bottom of the table. 
The problem of spuriosity of the Fermi matrix elements 
is here even more acute than in the case of Qvf3(i. In the 
absence of isospin violation, all Fermi matrix elements 
for 2j//9/3 should be exactly zero. Breaking of isospin is 
present in all calculations (ISM, QRPA, IBM-2, DFT, 
and HFB). Within the model space and in the closure 
approximation that we are using, we have a large break- 
ing for nuclei in which protons and neutrons occupy the 
same major shell and zero breaking in the others. The 
small value ~ 0.02 in Table IXVl is an indication of our nu- 
merical accuracy in calculating overlap of wave functions. 
From the dimensionless matrix elements in Table IXVl and 
the values of A we calculate the values of ImeC^M^^'^^l 



given in Table IXVI I In constructing this table we have 
taken into account only the GT matrix elements, since, 
as mentioned above, the IBM-2 F matrix elements are 
largely spurious in nuclei where protons and neutrons 
occupy the same major shell. 

We investigate two choices of Aqt- The first choice is 
that taken from Ref. ^l| or estimated by the systematics. 
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Act = 1.12^1/2 MeV, where A without tilde denotes the 
mass number. In cases where transitions between spin- 
orbit partners dominate, one expects the SSD approxi- 
mation to be appropriate. Our second choice is SSD for 
4oZr, 42M0, 46Pd, and 48Cd, where the dominant transi- 
tion is 59/2-57/21 E^iid 6oNd where the dominant transition 
is hii/2-hg/2. In the same table we also show the values 
of the matrix elements in the ISM without the closure 
approximation (38| . The ISM calculation are all in nu- 
clei in which protons and neutrons occupy the same ma- 
jor shell. By comparing these calculations with those in 
IBM-2 with the Fermi matrix elements set to zero we see 
that the two calculations have the same behavior with 
mass number but differ by a factor of approximately 2. 
The last columns in Table IXVII gives the values of the 
matrix elements lAf^/'^l extracted from experiment Q. 



□ from experimental T1/2 


(ISM) 


OgA 


= 1.25^*°'^ 


from experimental r 


,2 (IBM-2 


•gA 


= 1.25A-"'« 


□ H 








- & f=)\i_ 










\r** — 






Ca Ge Se 


Zr Mo Cd 


Te 


Xe t^d 



■40 60 80 100 120 140 160 

Mass number 



Figure 13. (Color online) Value of gA,eff extracted from ex- 
periment for IBM-2 and the ISM. 



If we write the matrix elements M21, as 

Af2t^= (^^)'m2., (38) 

where {gA,ef f I Qa) = (? is a quenching factor, by com- 
paring the experimental values M2u,exp with the calcu- 
lated values (or the experimental half-lives with those 
calculated using PSFs of jH) we can extract the values 
of gA.eff- These are given in Table IXVIII and Fig. [TH] 
for' IBM-2 (GT) and the ISM. As mentioned in Sec. HIl 
the renormalization of gA to gA.ef / is due to two main 
reasons: (1) limitation of the model space in which cal- 
culations are done and (2) omission of non-nucleonic de- 
grees of freedom (A,iV*. ...). As a result, one expects 
gA,eff to have a smooth behavior with A to which shell 
effects are superimposed. We see from Fig.[T3]that this is 
approximately the case if we assume SSD in 4oZr, 42 Mo, 
48Cd. and eoNd. This is consistent with previous analy- 
ses |85l.[86j. The smooth behavior can be parametrized as 
c,^BM-2 ^ 1 269^-'^, with 7 = 0.18 for IBM-2 (GT). This 
gives for the neutron (A=l) the free value. The same 
type of analysis can be done for the ISM. The values of 
gA,eff extracted by comparing the calculated and exper- 
imental matrix elements are also shown in Table IXVIII 
and Fig.[I31 We see that gA.eff in the ISM has the same 
behavior as in IBM-2, except for larger value. It can be 
parametrized as 5^^/ = 1.269^"''' with 7 = 0.12. In 
Ref. m the value 0.93 was used for ^SCa, ^^Qe and ^"^Se 
and 0.71 for "Oxe and "6Xe. 

The question of how to extract gA,eff in the QRPA 
has been the subject of many investigations [ll| . In this 
case gA,eff can be extracted either from 2^/3(3 or from 
single-/? decay jsgj. We do not discuss this extraction 
here but simply note that the values extracted are similar 
but larger than those in Table IXVlI] and Fig. [HI 



Values of {gA,effY can also be extracted from single-/? 
decay or electron capture using a Fermi-surface quasipar- 
ticle (FSQP) model [90j] where 

{gA,eff?=gffgff (39) 

is the product of gl^^ for the transition from even-even 
to odd-odd nuclei and g'j'^^ for the transition from odd- 
odd to even-even nuclei. The values obtained in this way 
(9H are also similar to those in Table IXVIII and Fig. [TS] 
Finally, very recently, values of 5^ ef'f\vf3p have been ex- 
tracted from a 21^/3/3 calculation without the closure ap- 
proximation for I28,i30xe_i.i28,i30xg decay with similar 
results [t^I. As one can see from the discussion in the 
paragraphs above, the extraction of the actual value of 
gA.eff is highly dependent on the model calculations and 
the assumptions made. All extractions, however, indi- 
cate values of gA^ef/^"^ the range 5^^^ ^ 0.57 — 0.90 
and 5^^y^ ^ 0.35 — 0.71 depending on mass number A 
and on the SSD or CA approximation, with decreasing 
trend with A. 

It is of considerable interest to analyze the impact that 
the quenching of gA to gA,eff observed in single-/? and 
2vl3p decay may have on Qvjifi. The question of whether 
or not the quenching of gA is the same in 2vj3(3 as in 
QvPP is the subject of debate, since only the states 1+ 
and 0^ in the intermediate odd-odd nucleus contribute 
to 2zy/3/3, while all multipoles contribute to Oi^/3/3. Two 
hues of thought have been considered : (1) Only GT (1+) 
is quenched and other multipoles are not. (2) All multi- 
poles are equally quenched. The experimental informa- 
tion on higher multipoles is meager, with only some hints 
coming from muon capture. The contribution of differ- 
ent intermediate states to Oj^/3/3 decay in ^°°Mo was 
investigated in Ref. jl3| within the framework of QRPA- 
Tii. It was found that the contribution of 1"*" is sizeble, 
of opposite sign of that of the other multipoles, and very 
much parameter (5pp) dependent. In view of this siz- 
able contribution, even if the other multipoles are not 
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Table XVI. Calculated values of 2!//3^ matrix elements in IBM-2 with qa = 1-269 and the ISM with qa = 1-25. 



A 



A(MeV) 



nieC M 



(21.) 





aCA 
A-GT 


ISSD 

Act 


IBM-2 




ISM^ 


exp^ 






CA 
GT 


SSD 
GT 




CA 

exp 


SSD 
exp 










48 


7.72^ 




0.10 




0.05 


0.038(3) 




76 


9.4l3 




0.24 




0.15 


0.118(5) 




82 


10.l3 




0.18 




0.15 


0.083(4) 




96 


ii.o" 


2.20 


0.10 


0.51 




0.080(4) 


0.075(4) 


100 


11.2 


1.69 


0.13 


0.89 




0.206(7) 


0.185(6) 


110 


11.8 


1.89 


0.13 


0.80 








116 


12.1 


1.88 


0.10 


0.63 




0.114(5) 


0.106(4) 


124 


12.5 




0.12 










128 


12.5^ 




0.15 






0.044(6) 




130 


13.33 




0.13 




0.07 


0.031(4) 




136 


i3.r 




0.11 




0.06 


0.0182(17) 




148 


13.6 




0.05 










150 


13.7 


1.88 


0.06 


0.42 




0.058(4) 


0.052(4) 


154 


13.9 




0.07 










160 


14.2 




0.11 










198 


15.8 




0.03 










Reference 


M- 














Reference 
















Reference 

















Table XVII. Value of QA.eff extracted from experiment. 



Nucleus 


ri/2,oxp(10'' yr)I 
exp 


Ti/2(10^*^ yr) 
IBM-2 




gA.eJJ 

IBM-2 




9A,eff 
ISM 






CA SSD 
GT GT 


CA 
GT 




SSD 
GT 





«Ca 


44+6 


2.30 




0.61(2) 




0.90(3) 


'"Ge 


1500 ± 100 


144 




0.71(1) 




0.90(2) 


«2Se 


92 ±7 


7.68 




0.68(1) 




0.74(2) 


^"Zr 


23 ±2 


5.31 


0.187 


0.88(2) 


0.38(1) 




"°Mo 


7.1 ±0.4 


6.46 


0.117 


1.24(2) 


0.46(1) 




"«Cd 


28 ±2 


14.5 


0.306 


1.08(2) 


0.41(1) 




128rpg 


1900000 ± 400000 


65600 


1170 


0.55(3) 








680li?° 


15.5 




0.49(2) 




0.67(3) 


"^Xe 


2110 ± 250^ 


23.0 




0.41(2) 




0.57(2) 




8.2±0.9~ 


3.21 


0.048 


1.00(3) 


0.35(1) 





Reference 
^ Reference 



quenched, there is going to be an effect coming from the 
l"*" multipole. 

In order to investigate the possible impact of quenching 
of gA , we present in Table IXVIIII the predicted half- lives 
under the assumption of "maximal quenching" in which 
all multipoles are quenched with 



^IBM-2 



(40) 



In the same table we show for comparison the un- 
quenched values with ~ 1.269 discussed in Sec.|IT] and 
shown in Table IXIIII We observe that while the ISM un- 



quenched values are a factor of approximately 2 smaller 
than IBM-2 values, the quenched values of both calcu- 
lations are similar, since the smaller calculated matrix 
elements in the ISM are compensated in part by a larger 
value of gA,eff- (This statement is not correct in ^^Ca 
where the simple parametrization 1.269 A"''' fails.) It ap- 
pears therefore that the difference of a factor of 2 in the 
calculated nuclear matrix elements in IBM-2 and the ISM 
is simply due to the difference in the size of the model 
space and thus in the renormalization of gA- 

Another important question in this context is whether 
or not gv is quenched. From the conserved vector current 
(CVC) we expect gv not to be quenched, at least as far 
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as the contribution of (A,iV*,...) is concerned. On the 
other side, the size of the model space certainly affects 
the Fermi matrix elements, through the overlap of the 
initial and final wave functions and their isospin purity. 
Thus, if one defines 

/ gv,eff,l3 \ _ \Mea:p^i3\ , , 

I 9v J ^ ' 

where gy = I, one may reasonably expect a quenching of 
the Fermi matrix elements as well. Whether or not the 
quenching factor for gv,eff,i3 is the same as for gA,eff.p 
is not clear. We are currently investigating this question 
within the context of IBM-2. Within this model, it ap- 
pears also that the question of isospin violation can be 
dealt with by means of a quenching of the Fermi matrix 
elements. The question of how to project into states of 
good isospin was investigated years ago by introducing 
the concept of F-spin (isospin of the pairs) jssl . p. 134]. 
Because of the complexity of the problem, we do not dis- 
cuss it here but defer it to a subsequent publication. In 
the columns "maximally quenched" in Table IXVIIII we 
have assumed equal quenching both for gy and gA and 
thus no quenching in the ratio gy/gA- This assump- 
tion introduces an additional error of about 10% in the 
quenched calculation. 

In conclusion. Table IXVIIII gives ranges of expected 
half-lives based on IBM-2 and the ISM calculations for 
unquenched, gA = 1.269, gv = 1, and "maximally 
quenched" values of gA.eff, gv,eff- The actual situa- 
tion may in fact be in between these two extreme values. 
Similar analyses have been done within the QRPA ex- 
cept that a quenched value gA,eff = 1-0 is used while 
gv.eff = 1 is unquenched fl3| . 



IV. CONCLUSIONS 

In this article we have presented a consistent evalua- 
tion of nuclear matrix elements in and Oi/hfi/S decay, 
Sec.|IT] and 2z^/3/3 decay, Sec. lIIII within the framework of 
IBM-2 in the closure approximation. All calculations can 
be done simultaneously by replacing the neutrino poten- 
tial v(j)) as summarized in Appendix A. While the clo- 
sure approximation is expected to be good for and 
QvhPP decay since the virtual neutrino momentum is of 
order 100 MeV/c and thus much larger than the scale of 
nuclear excitations, it is not expected to be good for 2vl3j3 
decay where the neutrino momentum is of order of few 
MeV/c and thus of the same scale of nuclear excitation. 
Furthermore, for 2i'/3/3, the single state dominance may 
be a better approximation. Hence the 21/(3/3 calculation 
in Sec. IIIII should be viewed only as an estimate. 

By using the Ov(3l3 matrix elements and phase-space 
factors of Ref. Q, we have calculated the expected Oiyf3l3 
half-lives in all nuclei of interest with gA = 1.269 and 
gv = l; given in Table IXlIIl and Fig. [51 This is the main 



result of this paper, and should be compared with other 
calculations, QRPA, ISM, DFT, and HFB, with the same 
(or similar) values of gA = 1.25 — 1.269 and gv ^ I- 

Finally, in Sec. IIIII we have examined the impact that 
a quenching of gA may have on Oi^/3/3 decay and reported 
in Table IXVIIII results of a quenched calculation with 
the quenching factor extracted from 21/(3/3 decay. This 
calculation is speculative since we have no experimental 
information to confirm whether or not quenching is the 
same for all multipoles in the intermediate nucleus. Our 
assessment is that, while in the unquenched case the sit- 
uation is such that current (GERDA and CUORE) and 
planned experiments may reach accuracies to detect at 
least the inverted hierarchy of Fig. [TUl in the quenched 
case only the degenerate case can be detected in the fore- 
seeable future. The results presented here point to the 
necessity of further studies and refinements, the crucial 
ones being: (i) an improved treatment of the Fermi ma- 
trix elements, (ii) an improved treatment of SRC, and 
(iii) the determination of the quenching factors gA,eff 
for all multipoles in 0i'/3/3 decay. The latter is of impor- 
tance not only for IBM-2 but also for all other model 
calculations. 
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VI. APPENDIX A: NEUTRINO POTENTIALS 
AND THEIR RADIAL INTEGRALS 

The neutrino potentials used in this article are given in 
Table IXiXl The function H{r) is the Fourier-Bessel trans- 
form of 2n~'^[p{p + A)]~^ , and is given in Appendix 
2, and in Eq. (19) of Q. It does not have an explicit 
form. We note, however that, when the closure energy A 
goes to zero, then v{p) = 2tt~^p~'^, and its Fourier-Bessel 
transform becomes the Coulomb potential, H{r) — > 1/r. 
The neutrino potential is a long-range potential, since 
the mass of the exchanged particle is very small. The 
situation is opposite in the case of heavy neutrino ex- 
change. In this case, the mass of the exchanged particle 
is very large and thus the potential is a contact interac- 
tion 5{r)/r'^ . For 2v/3(3 decay, the potential does not have 
a radial dependence and thus it is a contact interaction in 
momentum space. The values of A used in this article are 
given in Table IX VI I The radial integrals of the neutrino 
potential are best calculated in momentum space using 
the Horie method as discussed in the Appendix A of 
fl5| . with harmonic oscillator single-particle wave func- 
tions with oscillator parameter v = Muj/ti, where M is 
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Table XVIII. Predicted half-lives in OufiP decay with unquenched and maximally quenched qa, sjf^/^ and (/A,e// obtained 
from 2i//3/3 decay. 

Decay r°;,(10^Vr) 

IBM-2 ISM 



Decay 


unquenched 


maximally quenched 


unquenched 


maximally qi 


48Ca^48Ti 


1.03 


16.8 


13.9 


89.2 




1.45 


32.8 


8.65 


69.1 




0.52 


12.4 


2.22 


18.5 


^'^Zr^^'^Mo 


0.77 


20.5 






ioomo->iooRu 


0.46 


12.5 








1.60 


47.1 






ii6Cd^ii6Sn 


0.78 


24.0 








0.91 


29.2 


2.73 


27.6 


i28Te^i28xe 


8.53 


281 


33.5 


344.4 


iso^g^iaoxe 


0.44 


14.5 


1.70 


17.6 


"''Xe^i^^Ba 


0.62 


21.5 


2.39 


25.3 




2.54 


97.4 






isoNd^isogj^ 


0.30 


11.0 








5.34 


201 






160Q^^160i)y 


0.80 


31.0 






198p^^l98jj 


3.77 


170 







the nucleoli mass. In this article, we take ly = i^qA 
where A is the mass number and i/q = 0.994 fm^^. 



Table XIX. Neutrino potentials used in this article. 



Transition 


V{r) 


v{p) 




H{r) 


2 1 
vr p(p+A) 




1 S(r) 


2 1 

7r m^nip 


2ul3l3 


1 


S{p) 



VII. APPENDIX B: SINGLE-PARTICLE 
ENERGIES AND STRENGTH OF INTERACTION 

In order to calculate the pair structure constants 
we need the single-particle and single-hole energies and 
strength of interaction. We give in Tables IXXllXXIII the 
single-particle and single-hole energies used in this arti- 
cle. We generate the pair structure constants by diag- 
onalizing the surface delta interaction (SDI) in the two 
identical particle states, pp and nn. The strength of the 
(isovector) interaction, Ai, is also given in Tables IXXI 
IXXIII It is obtained by fitting the 2+ 0+ energy differ- 
ence in nuclei with either two protons (proton holes) or 
two neutrons (neutron holes). For ^^Ca— "^^Ti decay, we 
need also the strength of the interaction in the I/7/2 shell, 
given by Ai = 0.510 MeV. The calculation of the pair 
structure constants can be improved by a better choice 
of the interaction and of the single-particle energies. We 
have tried different choices of the single-particle energies 



Table XX. SDI strength values A\ and single-particle and 
single-hole energies (in MeV) in the N,Z — 28 — 50 shell. 
The energies are taken from the spectra of ^^Cu for proton 
particles, from isotones A*' = 50 for proton holes, and from 
the spectra of ^^Ni for neutron holes. 





Protons 


Protons 


Neutrons 


Orbital 


(particles) 


(holes) 


(holes) 




Ai = 0.366 


Ai = 0.264 


Ai = 0.280 


2pi/2 


1.106 


0.931 


1.896 


2p3/2 


0.000 


2.198 


3.009 


l.f5/2 


1.028 


2.684 


2.240 


139/2 


3.009 


0.000 


0.000 



Table XXI. SDI strength values Ai and single-particle and 
single-hole energies (in MeV) in the N, Z = 50 — 82 shell. 
The energies are taken from the spectra of ^^^Sb for protons 
particles, from the spectra of ^''^Tl for proton holes, from the 
spectra of '^Zr for neutron particles, and from the spectra of 
^''^Sn for neutron holes. 





Protons 


Protons 


Neutrons 


Neutrons 


Orbital 


(particles) 


(holes) 


(particles) 


(holes) 




Ai = 0.221 


Ai = 0.200 


Ai = 0.269 


Ai = 0.163 


3Sl/2 


2.990 


0.000 


1.205 


0.332 


2rf3/2 


2.690 


0.350 


2.042 


0.000 


2rf5/2 


0.960 


1.670 


0.000 


1.655 


I37/2 


0.000 


2.700 


2.200 


2.434 


1^11/2 


2.760 


1.340 


2.170 


0.070 



and included the variation of the corresponding radial 
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Table XXII. SDI strength values Ai and single-particle ener- 
gies (in MeV) in the A'^ = 82 — 126 shell. The energies are 
taken from |93 | for neutron particles and from the spectra of 
^°*Pb for neutron holes. 





Neutrons 


Neutrons 


Orbital 


(particles) 


(holes) 




Ai =0.147 


Ai =0.150 


3pi/2 


2.250 


0.000 


3P3/2 


1.500 


0.900 


2/5/2 


2.600 


0.570 


2/7/2 


0.000 


2.340 


l/l9/2 


2.450 


3.410 


1*13/2 


2.800 


1.630 



VIII. APPENDIX C: PARAMETERS OF THE 
IBM-2 HAMILTONIAN 

A detailed description of the IBM-2 Hamiltonian is 
given in (IJ and (U. For most nuclei, the Hamilto- 



integrals in the estimate of the sensitivity to parameter 
changes. 



nian parameters are taken from the literature (95l4l07 |. 
The values of the Hamiltonian parameters, as well as 
the references from which they were taken, are given in 
Table IXXIIII The quality of the description can be seen 
from these references and ranges from very good to excel- 
lent (see Fig. 3]). The only nuclei for which we have done 
new calculations are -^^Ti, ^^Zr, i24Sn, ^^^Xe , i^Qd, and 
i60j3y ^jjg jjg.^ calculations are done using the program 
NPBOS (H adapted by J. Kotila. They include ener- 
gies, B{E2) values, quadrupole moments, B{M1) values, 
magnetic moments, etc. The c alcul ations for ^^°Gd and 
160 have just been publi shed |l08| . A paper with those 
for *^Ti is in preparation |l09l |. The quality of these, as 
well as of the unpublished results for ^^Zr is equal to that 
of the results obtained previously (95l4l07| . For the semi- 
magic nuclei ^^^"^^''Sn and ^'^^Xe, we have obtained the 
parameters by a fit to the energy o f the low-lying states 
using the same procedure as in Ref. |103| for ^^^Sn, while 
^®Ca has been taken as doubly magic. This procedure is 
compatible with the generalized seniority (GS) scheme, 
which appears to be good for semi-magic nuclei, as exten- 
sively discussed in the 1980s for pairing plus quadrupole 
intera ction s, and as shown recently for realistic interac- 
tions fiTot . 
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Table XXIII. Hamiltonian parameters employed in the IBM-2 calculation of the final wave functions along with their references. 

Tvr i f c f (0) (2) (4) (0) (2) (4) 

■'^^Ti,^ 1.11 1.11 -0.20 -0.30 -0.70 1.00 1.00 1.00 

■^''Ge [95] 1.20 1.20 -0.21 1.00 -1.20 -0.05 0.10 -0.05 

'■'^Se [96| 0.96 0.96 -0.16 0.50 -0.90 -0.10 

^^Se [96| 1.00 1.00 -0.28 1.14 -0.90 -0.10 

*^Kr [97| 1.15 1.15 -0.19 0.93 -1.13 -0.10 -0.10 

^^Zr^ 1.00 1.00 -0.20 -2.20 0.65 0.17 0.17 0.33 

^^MojgSl 0.73 1.10 -0.09 -1.20 0.40 -0.10 0.10 -0.10 -0.50 0.10 

""Mo [98| 0.55 1.00 -0.06 -1.20 0.40 -0.10 0.10 -0.10 -0.60 0.20 0.10 

"°Ru [99] 0.89 0.89 -0.18 -1.00 0.40 0.60 0.09 -0.13 

""Pd [100| 0.78 0.60 -0.13 0.00 -0.30 0.20 0.04 0.00 -0.26 -0.29 -0.30 -0.26 -0.29 -0.03 

"°Cd [lOlJ 0.92 0.92 -0.15 -1.10 -0.80 1.10 0.109 1.10 0.07 -0.17 0.16 

"''Cd [102| 0.85 0.85 -0.27 -0.58 0.00 -0.18 0.24 -0.18 -0.15 -0.06 

""^Sn [103| 1.32 -0.50 -0.22 -0.07 -0.06 0.04 

^^''Sn^ 1.10 -0.30 -0.16 -0.20 0.30 0.02 

^2*Te [i02| 0.82 0.82 -0.15 0.00 -1.20 -0.18 0.24 -0.18 0.10 

i^^Te [102| 0.93 0.93 -0.17 0.50 -1.20 -0.18 0.24 -0.18 0.30 0.22 

i2**Xe [104| 0.70 0.70 -0.17 0.33 -0.80 -0.18 0.24 -0.18 0.30 

^^°Te [102| 1.05 1.05 -0.20 0.90 -1.20 -0.18 0.24 -0.18 0.30 0.22 

"°Xe [104| 0.76 0.76 -0.19 0.50 -0.80 -0.18 0.24 -0.18 0.30 0.22 

^Xe*" 1.31 -0.04 0.01 -0.02 
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"®Ba [104| 1.03 1.03 -0.23 1.00 -0.90 -0.18 0.24 -0.18 0.30 0.10 

^***Nd [105| 0.70 0.70 -0.10 -0.80 -1.20 -0.12 0.24 0.90 0.40 0.20 

"*Sm [105J 0.95 0.95 -0.12 0.00 -1.30 -0.12 0.24 0.90 0.05 

^^°Nd [105| 0.47 0.47 -0.07 -1.00 -1.20 -0.12 0.24 0.90 0.40 0.20 

^^°Sm [105| 0.70 0.70 -0.08 -0.80 -1.30 -0.12 0.24 0.90 0.05 

^''''Sm [105| 0.43 0.43 -0.08 -1.10 -1.30 -0.12 0.24 0.90 0.05 

^^*Gd [105| 0.55 0.55 -0.08 -1.00 -1.00 -0.12 0.24 0.90 -0.20 -0.10 

"°Gd [108| 0.42 0.42 -0.05 -0.80 -1.00 0.08 0.08 0.08 -0.20 -0.10 

"°Dy [108| 0.44 0.44 -0.06 -0.80 -0.90 0.08 0.08 0.08 -0.05 -0.15 

i^^^Pt [106| 0.58 0.58 -0.18 1.05 -0.80 -0.10 0.08 -0.10 0.00 0.02 0.00 

"*Hg [107| 0.55 0.55 -0.21 1.00 -0.40 0.08 0.37 0.25 0.16 

Parameters fitted to reproduce the spectroscopic data of the low-lying energy states. 

^ GS parameters fitted to reproduce the spectroscopic data of the low-lying energy states. 
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